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1. Introduction 

The notions of 'free probability' and 'freeness' of non-commutative random variables 
were introduced by Voiculescu in the study of certain algebras of bounded operators 
[25] . The word freeness in this context means a kind of statistical independence of 
operators. The algebraic concept of freeness of random variables has a natural realization 
by random matrices in the limit of infinite matrix dimension [23] ; this realization is what 
we study, develop, and apply in the present paper. 

A central tool of the free probability formalism is the so-called /^-transform, which 
resembles the logarithm of the characteristic function for commutative random variables. 
Voiculescu [23] defined it by the functional inverse of the average trace of the resolvent 
operator. A second approach to the subject is due to Speicher [T9|, who expressed 
the moments of the random matrix directly in terms of the Taylor coefficients of the 
i?-transform. Speicher's concept of non-crossing partition is a mathematical expression 
of the dominance of planar Feynman graphs (using physics parlance) in the large- 
limit. In the present paper we will encounter both approaches: the analytical one of 
Voiculescu, and the combinatorial one of Speicher. 

Our long-term goal is a comprehensive description of the spectral correlation 
functions and, ultimately, a proof of the universality hypothesis which is expected for 
certain random matrix ensembles in the large- A^ limit. Although the i?-transform is a 
powerful tool to tackle density-of-states questions, it is fair to say that free probability 
theory has not yet contributed much to our understanding of the universality of spectral 
correlation functions at the scale of the level spacing. 

Bearing this in mind, we now change subject and turn to the so-called 
supersymmetry method, by which we mean the technique of integration over commuting 
and anti-commuting variables pioneered by Wegner [26] and Efetov [10]. In its original 
formulation (using the Hubbard-Stratonovich transformation) this method was limited 
to Gaussian disorder distributions. Nonetheless, with this limitation it has enjoyed 
great success in producing non-trivial results for a number of physics problems, e.g., 
level statistics of small metallic grains, localization in thick disordered wires, scaling 
exponents at the Anderson transition, etc. In the present paper we will take advantage 
of a recent variant, called superbosonization, which makes it possible in principle to 
treat a class of disorder distributions much wider than the Gaussian one. 

Since their inception in the 1980 's, free probability theory and the supersymmetry 
method have coexisted with little or no mutual interaction. Forecast in a prescient 
remark by P. Zinn- Justin [30], the message of the present paper is that a new quality 
emerges when the two formalisms are combined. More specifically, we will show that the 
characteristic function of the probability law of the random matrix ensemble - an object 
of central importance to superbosonization - has a large- A^ limit which is determined 
by the /^-transform. This result paves the way for a number of applications. As a 
first application, we will illustrate the new method by demonstrating universality for a 
random matrix model of stochastic scattering. 
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1.1. Summary of results 

Our results make reference to the i?-transform, which we now introduce in more detaiL 
Consider the average trace of the resolvent operator, g{z) := limAr_>oo A^~^(Tr {z — H)~^), 
2; G C \ M , of the random matrix H. Voiculescu inverts the function z ^ g{z) to define 
the /^-transform as the regular part of g{z) = k ^ + R{k) = z . It has a power 
series R{k) = J2n>i whose coefficients c„ are called free cumulants. These are 

analogs (in the non- commutative setting) of the usual cumulants in that they are linear 
with respect to free convolution. 

In the present paper, we consider Utv- invariant probability measures with a 
density of the form 

:=e-^T''^(^)rfif , (1.1) 

where dH denotes Lebesgue measure on the linear space of Hermitian N x N matrices, 
and R 9 X I— )■ V{x) is analytic and uniformly convex. The focus of our analysis is the 
characteristic function 

Q{K) = j e^'^'^'dfiNiH) . (1.2) 

Note that Q{K) is invariant under conjugation K i— )■ g^^Kg hj g E GLn . 

Motivated by the supersymmetry method as reviewed in Section [2l we make a study 
of fl{NK) for linear operators K of the form 

p q ^ 

K = ^(Pa(S> ^a + '^A'^^b : (1-3) 
a=l b=l 

with p,q kept fixed in the limit — )■ oo. Here ipa,'^b ^ are vectors, while 
'^a ,'^b & (C^)* are co- vectors. The components of (pa are complex commuting variables 
or generators of the symmetric algebra S(C^), whereas the components of ipb are 
complex ant i- commuting variables or generators of the exterior algebra A(C^). The 
same statement, albeit with replaced by (C^)*, applies to ipa and ipb respectively. 
All inner products {if a, fb) for a,b = 1, . . . ,p are kept fixed in the large- limit. 
For such K we argue that the following holds: 

oo 

lim A^-^ In Q(NK) = V — Tr JsT" . (1.4) 

n=l 

We note that this formula makes perfect sense as long as the power series of the R- 
transform has infinite radius of convergence. The latter property is ensured by our 
assumptions on V. 

1.1.1. Related mathematical work. To put the result ( 11.41) into context, we now mention 
related mathematical work on the large- asymptotics of the following spherical integral 
(known in the physics literature as the Itzykson-Zuber integral): 

IxAk)= I e^^^(^-^"^-^)rf^, (1.5) 
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for the case of a rank-one projector 11. Let the eigenvalues Xi^n, ■ ■ ■ jX^^n of the 
Hermitian matrices be confined to a finite interval [a, b] and assume that the 
empirical measure 6{x — xj^n) converges weakly to a measure with support 

in [a, b]. Under these conditions the following is known. 

Collins [6] differentiates the scaled logarithm of the spherical integral n times at 
zero to show that 

hm iV-i— In IxAm ={n- 1)! , (1.6) 

i.e., he establishes convergence to the n^^ free cumulant (times a factorial). A stronger 
version of this result, 

„k oo 

lim N^^ In Ix^{Nk) = / R{k') dfc' = V — fc" , (1.7) 

TV^oo Jo ^ n 

was proved by Guionnet and Maida [15] under the condition that A; G C is small enough. 
(Notice that (11. 7p implies (II. 6p .) For k real and large, however, the authors of |15] obtain 
a different behavior, separated from the small-fc regime by a phase transition. 

These results have a bearing on (ll.4p for p = 1, q = because the integral (II. 2p can 
be done in two steps: fixing some set of eigenvalues for H we first do the integral over 
Utv orbits - that's precisely the spherical integral (II. 5p - and afterwards we take the 
average over the fluctuating eigenvalues. While it may seem puzzling at first sight that 
the authors of [15] find a phase transition whereas we do not. Section explains that 
there is no contradiction here. The assumption of eigenvalues strictly confined to an 
interval [a, b] would mean in our context that the confining potential V{x) is infinitely 
high outside of [a,b]. In contradistinction, we assume that V is both uniformly convex 
and analytic. In the latter setting the existence of a phase transition is ruled out on 
physical and mathematical grounds. We in fact argue that the limit in (II. 7p is an entire 
function of G C. 

1.1.2. Towards applications. By the general principles of invariant theory, the charac- 
teristic function fl{K) = fl{g~^Kg) {g G GLtv) lifts to a function Q{Q) of the super- 
matrix of GLTv-invariants 

\ / a,a'=l,...,p; b,b'=l,...,g 

The angular brackets still mean contraction of the vector with the co-vector. By the 
relation Tr fT" = STr (where STr denotes the supertrace) it follows from (ll.4p that 

oo 

lim N-^ In Q(NQ) = V — STr Q" . (1.9) 

/V— iro ' ^ n 



N^oo ^ — ^ n 

n=l 



The lifted function is input required by the superbosonization method (Section |2]); 
in the past this input was the link missing for applications. With (II. 9p established, 
we now have at our disposal a powerful new method for the treatment of random 
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matrix problems. In the present paper we illustrate the new method by demonstrating 
universality for a random matrix model of stochastic scattering. In particular, we point 
out that the pertinent large-A^ saddle-point equation for Q is Voiculescu's equation 
+ R{k) = z generalized to the (super-)matrix case: 

+ R{Q) = z ldj,\g . (1.10) 

While the present paper addresses only the case of unitary symmetry, our treatment 
is robust and readily extends to ensembles with orthogonal or symplectic symmetry. 

1.2. Outline 

An outline of the contents of the paper is as follows. Section |2] provides background 
and motivation by introducing the characteristic function Q{K) as a key object of the 
supersymmetry method. For the special case of K = ^ If = kU with A; G M and 
a rank-one projector 11, the large- asymptotics of Q{NK) is computed in Section [3] 
using Dyson Coulomb gas methods. Particular attention is paid to the fact that the 
asymptotics for small and large k match perfectly to give an answer which is smooth as 
a function of k. The fermionic analog K = ip ® with ant i- commuting ?/^, tp is treated 
by drawing on information from representation theory in Section HI By using standard 
perturbation theory in the large- limit we then develop in Section [5] a combinatorial 
description of the full superfunction VL{NK). The resulting formalism is applied to a 
model of stochastic scattering in Section [61 An outlook is given in Section [71 

2. Review: supersymmetry method 

We begin the paper with a concise review of the supersymmetry method and, in 
particular, of superbosonization. In this way, we shall introduce and motivate the Fourier 
transform VL{K) = J e^'^^^dfii\f{H) with K given by (11.31) . which is a superfunction with 
symmetries and the key object to be analyzed in the sequel. 

2.1. First steps 

Consider the Hermitian vector space with its standard Hermitian scalar product 
h : X — )■ C , which determines a C-antilinear bijection 

t:C^^(C^)*, v^v^ ■.= h{vr) , (2.1) 

between and its dual vector space (C^)*. Following standard physics conventions, 
we denote by the same symbol (dagger) the operation L i— )■ of taking the Hermitian 
conjugate of a linear operator L G End(C^). 

In the following the Hamiltonian H will always be a random Hermitian operator: 

End(C^) 3 H = H\ (2.2) 

distributed according to some probability measure fij\[ . Our goal is to study the spectral 
correlation functions which are defined as averages with respect to fiN ■ For this purpose 
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we consider the characteristic polynomial z i->- Det(2; — H) associated with H. The 
supersymmetry method allows us to compute expectations of products of ratios of such 
polynomials, and hence of products of resolvent traces Ha'^^ i^a '~ as follows. 

Let us denote the canonical pairing (C^)* ® ^ C (i.e., evaluation of a linear 
form if e (C^)* on a vector tp e C''^) by 

(fi ^ (p {(fi,(p) . (2.3) 

With the resolvent operator (z — H)~^ for z G C\]R we associate a holomorphic function 
7 : (C^)* X ^ C by 

7(^, (f) = e-<^''^^-^'^> . (2.4) 

Now let C (C^)* X be the graph of the K-hnear mapping 

C^^(C^)*, v'^-iV, s := sign(am2) G {±1} . (2.5) 

Thus is the real vector space of all pairs {(p,(p) = {—isip\(p) for ip G C^. The 
Gaussian 7 decreases rapidly along m . Indeed, 

d\t{lfi,ipz- H(p) ^ Dlth((p,-is(z - H)(p) ^ \Jm z\ h((p,(p) >0, (2.6) 
so we may integrate 7 along . By a standard formula for Gaussian integrals we have 

^-(^,^z-H^} ^ i3et-i(^ - H) , (2.7) 

where the integral is over m with (iR-valued) Lebesgue measure normalized by the 
condition Jy^ Q-i'f^f) — 1. (This measure is not made explicit in our notation.) 
Expressing Tr (2; — H)"^ as a logarithmic derivative, 

we see that we need a Gaussian integration formula for Det(2; — H) (where 2; G C) in 
addition to that for reciprocals Def~^(2; — H). Such a formula can be had by replacing 
commuting variables </? by anticommuting variables ^0 ; i.e., we view 

{i^,^^^z-H^^;)e^'{C'' ^iC")*) (2.9) 

as a quadratic element of the exterior algebra generated by the direct sum © (C^)*. 
The precise meaning is this. Let {e^} be a basis of and {e*} be the dual basis of 
(C^)*. LetJ : © (C^)* ^ A(C^ © (C^)*) be the canonical embedding; or simply 

put, view ipi = i{ei) and = ^(e*) as anticommuting variables or generators of the 
exterior algebras A(C^) and A((C^)*) respectively. Then we define {i/jjtjjz — Ht/j) to 
be the element of A^{C^ © (C^)*) given by 

{^,^z- H^) V. . V^, {e\ {z - H) sj)^^ . (2.10) 

By exponentiating this expression we get a Gaussian element in (the even part of) 
the full exterior algebra: 



G 



(f) A''=(C^ffi(C^)*) . (2.11) 
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The Berezin integral f ^ J f for / G A(C''^© (C'^)*) is, by definition, projection on the 
one-dimensional subspace A^^(C^ © (C^)*) of top degree. In the case of our Gaussian 
integrand the result of this projection is known to be proportional to the determinant of 
the operator z — H. We normalize the Berezin integral in such a way that the constant 
of proportionality is unity: 

^{i,,^z-Hi,) ^ _ _ ^2.12) 

To summarize the above, we have two Gaussian integration formulas: Eq. (12.12^ for the 
secular determinant Det(z — H), and Eq. (12. 7p for its reciprocal. 

By multiplying these formulas, averaging the result with the given probability 
density dfiN , and interchanging the order of integrations, we obtain 

Det{w,-H) ^^^^^^ ^ I ^](^©^ + V;©^)e-""'<^■^>+"'^<'^•'^>(2.13) 



where i^i G C , iuq G C \ M , and Vt is the characteristic function 

n{K) = j e'^'^'^'di^NiH) . (2.14) 

Note that the formula (12.131) requires the evaluation of Vt for 

((^Vi + V'Vi)ei (e^ ■) , (2.15) 

where ^9* : — t- C and (pj : (C^)* C are the linear coordinates associated with the 
bases {cj} of and {e-'} of (C^)*. Thus the operator K is an endomorphism of 
with coefficients in the tensor product 

S(C^©(C^)*)©A(C^©(C^)*) (2.16) 

of the symmetric and exterior algebras of © (C^)*. It will be important that the 
numerical part oiK = Lp®ip + ip®ip has finite rank. 

The relation (12.131) transfers the integral over N x N random matrices H to an 
integral over the variables ip constituting the bilinear K. This transfer will be a step 
forward if we can calculate the function Q{K) or, at least, gather enough information 
about it. For the case of a Gaussian probability measure /iat, the Fourier-Laplace 
transform Q{K) is also Gaussian. The supersymmetry formalism then takes its course 
and delivers results quickly. However, using the traditional version of the supersymmetry 
method one did not know how to proceed in the general case of non-Gaussian /i^v • 



2.2. Superbosonization 

One way to proceed, as we shall now review, is to make a symmetry assumption about 
Q{K). Let V := (C^)* © C^. Then e^"^ for fsT = © ^ + © is a superfunction 
/ : K — 7- /\{y*), and by making the identifications ip"^ = dip^ and ipi = dipi we may 
regard / = e"'"''^^ as a holomorphic differential form on the complex vector space V. 
Let now //at be invariant under conjugation H gHg~^ by the elements g of some 
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compact group G. Then by integrating / against we obtain a differential form 
Q : V ^ A(y*) which is G-equivariant, i.e., 

Q{v) = gn{gv) , (2.17) 

where {g,v) (-)■ gv and fi) i— )• gQ are the natural G-actions on V resp. A{V*). Later 
we will write this equivariance property more intuitively as Q{K) = Q{gKg~^). 

In such a setting, the superbosonization method offers a reduction step which is 
available [H] for the classical Lie groups G = \Jn , On , and USp^ . For each of these 
groups, the algebra of G-equivariant differential forms on V is generated by (the dual 
of) the Z2-graded vector space W = Wq © Wi of quadratic G-invariants [16j. After 
lifting the differential form Q : V ^ A(y*) to a superfunction Q : Wq ^ A(W^j*), the 
step of 'superbosonization' transfers the integral on the r.h.s. of (12.131) to an integral of 
the lifted superfunction over a (low-dimensional) Riemannian symmetric superspace. 

We now reproduce from [TT] the details for the case of G = Un , with a few 
notational adjustments to fit the present situation. By immediate generalization of 
(I2.13P we have 

Ut=i^et{w,,, - H) 
nL.Det(^o,.-i/) '^"^""^ - 

/ ^(5^^<^a(^ar) + 5Z^^6(^6,-))e-^"'«-<^»'^»>+^"'^'''<'^-'^''>, (2.18) 

where the (/^-integral is over the real subspace 

= -iSaV^l , = sign(Jmwo,a) , a = l,...,p. (2.19) 

While intending to specialize to the case of p = q later, we here describe the general 
case p ^ q for a clear exposition of the formalism. 

To simplify the notation, it is convenient to regard the vectors (pi,...,(pp as the 
components of a linear mapping ip from to : 

^ := (y^i, . . . , </?p) e Hom(CP, C^) . (2.20) 

Similarly, we view (p := {(fi, . . . , (fp) as a linear mapping 

^ G Hom(C^, C*') . (2.21) 

Using the same conventions on the anti-commuting side, we write our integral as 
nLiDetK,,-.ff) 
riL.DetK.-i/) """^"l 

+ ^tfj) e-^' iv^o-fi+^^i^) _ (2.22) 



Here Tr = Tr^iv , and wq = diag(wo,i, . . . , wo,p), wi = diag(wi,i, . . . , wi,q) are diagonal 
operators. The integral is over 

^=-is(^t^ s = diag(si,...,Sp) . (2.23) 
It is evident that the integrand on the r.h.s. of (12.221) has the invariance property 

f{ip,ip,ip,^) = f{gip,^g-^,gip,^g-^) (2.24) 
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for g ^ Vn and hence, by holomorphic continuation, for g G GL^v • This imphes that 

X a 
T y 



there exists [.17] a (hfted) function f{Q) of a supermatrix Q = I ) such that 



/I f f* |=/(^.^,^.,^.). (2.25) 

Precisely speaking, / : Wq — )■ /\{Wi) is a function on the Z2-even vector space 
Wq = End(C^) © End(C^) (the diagonal blocks of Q) with values in the exterior algebra 
of the dual of the Zg-odd vector space Wi = Hom(CP, C«) © Hom(C^ C^) {the off- 
diagonal blocks of Q). The lift f{Q) turns into the given function f{ip,ip,ip,ip) upon 
substituting for Q the quadratic GLAr-invariants, 

Q=( ' '']^( . (2.26) 

\T y j yipip i>i> J 

The ensuing step of superbosonization exploits the GLTv-symmetry of the integrand 
to implement a reduction: it transfers the integral over the variables if, ip to an integral 
over the supermatrices Q : 

f{^,^,i:,i^) = c,,,j SDet^(g) /(g) . (2.27) 

The integral on the right-hand side is over 

H; X Ug C End(C'') X End(C^) , (2.28) 

where Uq = U(C'') is the unitary group of and 

:= {-IsM \ M = >0} (2.29) 

is a space isomorphic to the positive Hermitian p x p matrices (replacing the quadratic 
Utv" invariant v'V)- With a natural choice [17] of normalization for the Berezin 
integration form DQ the normalization factor is Cp^g = vol(U„)/vol(U„_p+g). In the 
important special case of p = q the Berezin integration form DQ is simply the product 
of differentials times the product of derivatives: 

Dg OC J] dXaa' n ^^ybb' H a a ' (2-30) 
a, a' b, b' a, b 

Finally, it should be stressed that the formula ( I2.27P is valid if and only if > p. 

Application of the superbosonization formula fl2.27p to Eq. fl2.22p yields the identity 

nLiDet(t.l,>-g) 



j DQ SDet^(g) n{Q) e'"^^' ""^ , (2.31) 

where STr wg = TT(Cp{woQ)—TTcq{wiQ). The function ^{Q) is a lift of the characteristic 
function Q{K) for K = ipip + ipip. 

In view of the result (12.310 our short-term goals should now be well motivated: 
the key object to understand is the lifted characteristic (super-)function VL{Q). If we 
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can control Q , results for the level correlation functions will follow from a large- 
asymptotic saddle analysis of the Q-integral. 

3. Coulomb gas argument 

Based on what is called the Dyson Coulomb gas, we are now going to study Q{K) for the 
rank-one case K = kU with G M and 11 the projector on a one-dimensional subspace 
of C^. A related situation has been investigated in the work of P. Zinn- Justin [291 130] . 
Collins [6], and Guionnet & Maida [I5]; we will comment on the literature as we go 
along. We begin by reviewing some basic material. 

3.1. Voiculescu R-transform 

As before, let /ijv be a probability measure (with or without invariance properties) for 
Hermitian N x N matrices H. Let 



:= J Tr (H^) rf/i^(if) (3.1) 

denote the n}^ moment of fi^ ■ We shall assume that m„^jv has a limit, 

m„ := lim m„ n , (3.2) 
and consider a generating function z ^ g{z) for these moments at A^ = cxd : 

oo 

g{z):=Y,m^z—' . (3.3) 

n=0 

This series (when it converges), or rather its analytic continuation 

g{z)= lim N-^ [ {z - H)-^d^iN{H) (z G C \ R) , (3.4) 

Af— >oo J 

is the Cauchy transform 

q(z) = / , du(x) = 7r~^ lim 3mq(x — ie) dx , (3.5) 

JrZ-X ' e^0+ 

of the large- A^ limit of the so-called density of states u. Let us now assume that u has 
compact support. Then there exists a number r > such that for z G C with |z| > r 
the power series (13. 3p converges, the derivative g'{z) does not vanish, and hence by 
the implicit function theorem the function z i-> g{z) has a local inverse. Because the 
expansion of g{z) around z = oo begins as A; = g{z) = z~^ + . . . , the power series for the 
inverse function will likewise begin as z = plus corrections. These considerations 
lead to Voiculescu's definition of the i?-transform ^23j : 

g{z) = k^k~^ + R{k) = z , R{k) = V c„ . (3.6) 

Thus the /^-transform k t— )■ R{k) is the function inverse to z i— ?■ g{z) with the pole k~^ 
subtracted. Under the assumption of compact support for v the power series for R{k) 
converges for sufficiently small k . The coefficients c„ are called free cumulants. 
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Freeness of random variables is defined in an algebraic way which will not 
be reviewed here. Suffice it to say that two random matrices A, B are free (in the 
limit = oo) if the probability law oi A + B remains unchanged under conjugation 
B I— )■ WB U hj any U G Uqo • Voiculescu proved [23] that the i?-transform is linear 
for free convolution, i.e., if A, B are free [with i?-transforms RA{k) resp. Rsik)], then 
RA+B{k) = RA{k) + RB{k). This linearity property parallels the fact that the cumulants 
of a sum of commutative random variables equal the sum of the cumulants, and it leads 
to the expectation that the i?-transform is closely related to the logarithm of the Fourier 
transform of fi^ , N ^ oo. 

Let us mention in passing that, since g{z) is also known as the Green's function, 
physicists by Zee's fancy [28] sometimes call b{k) := k~^ + R{k) the Blue's function. 

3.1.1. Examples. The Gaussian measure fi^ with density dfiN{H) oc e"^'^^^^ dH is 
called the Gaussian Unitary Ensemble (GUE) with C2 = 1. For this measure one has 
R{k) = k and solving the equation z = k^^ + R{k) = k~^ + for /c = g{z) one finds 

g{z) = l{z±V^^), (3.7) 

which gives Wigner's semicircle law di'{x) = (2tt)'^^\/A — x"^ dx . 

Another example (taken from recent work [18j by Lueck, Sommers and one of the 
authors, on the energy correlations of a random matrix model for disordered bosons) is 
this. Consider 

^(^) = Y^ = E^""- (3-8) 

n=l 

Thus all odd free cumulants vanish and the even free cumulants are all equal to unity. 
Solving Voiculescu's equation k~^ + R{k) = z for k = g{z) one obtains 

The density of states of this example has compact support but is unbounded due to an 
inverse cube root singularity at x = . 

3.2. Eigenvalue reduction of Q 

We now specialize to the case of a UAr-invariant probability measure for Hermitian 
operators H with density dfi]\f{H) = e~^^^^^^^dH . The goal here is to relate the R- 
transform R{k) to the large-iV limit of the characteristic function (11. 2p for K = kU 
with n a rank-one projector. Our approach will be similar to that of P. Zinn- Justin [30] 
based on the Harish-Chandra-Itzykson-Zuber integral. 

We start by diagonalizing the Hamiltonian if by a unitary transformation: 

H = g~^Xg, X = diag(xi, . . . ,XAr) , gel^N- (3.10) 
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Recalling that the Jacobian J{X) associated with this transformation is the square of 
the Vandermonde determinant: 

JiX) = llix.-x,)\ (3.11) 

i<j 

we cast the expression for the characteristic function in the form 

Q{K)=Cj,[ ([ eT^(^^^^-')f/(7')e-^T^^Wj(X)rf^a:, (3.12) 

where dg is a Haar measure for Uat . The normalization constant Cat is determined by 
the condition Q{0) = 1. 

Now let K = kU where A; G M and 11 is the orthogonal projector on some (fixed) 
complex line in C^. We are then faced with the inner integral 

The integrand depends on G Uat only through the projector 11 conjugated by g, 
and the set of all these projectors gllg~^ is in bijection with the projective space 
([^pN-i ^ Uiv/(Ui X Ujv_i) of complex lines in C^. By parametrizing gllg~^ in the 
eigenbasis of H as {gllg~^)ij = uiUj with Yl!j=i I'^iP — reduce our integral to 

1 ^k^(xgug~^)^g^ f e'^U-^M'du, (3.14) 

Jvn 

where du is a UTv-invariant measure for CP^^^. 

Now CP^^^ is a Kahler manifold with UAr-invariant Riemannian geometry, and 
the function 

/i : CP^"^ ^ Lie Utv , ^ ■ (Ui x U^v-i) ^ i^n^~\ (3.15) 

is a momentum mapping We observe that the expression kTT{XgIlg~^) in the 
exponent of our integrand is obtained by contracting fi with the Lie algebra element 
—ikX G Lie U^v • It follows that the integral is governed by the Duistermaat-Heckman 
localization principle pj. In other words, the integral can be computed exactly by 
performing the stationary-phase approximation (including the Gaussian fluctuations) 
for each of its critical points and summing the contributions. 

There are critical points; these are the points where gllg~^ is diagonal (with 
one diagonal element equal to unity and all others equal to zero). By computing the 
contribution from each point and taking the sum, we get 

r ^ kx 

/ e'^U^^\'^^^du = C,,^k-^^-'^Y.T^^ V ' (3-16) 

with some A^-dependent constant Ci^n ■ 
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3.3. Coulomb gas 

Based on the exact expressions ( ]3.12p -( l3.16p . our goal is to compute the large- 
asymptotics of A^~^ lnQ{NkIl). We begin by recalling |8j that, governed by 

J](x,-x,)2j]e-^^(^')da;z, (3.17) 

i<j I 

the eigenvalues Xi, . . . ,X]y distribute for — )• oo according to the equilibrium measure, 
z/, which is determined by minimizing Dyson's Coulomb gas energy functional: 

j V{x) dv{x) jj In |a; - y| dv{x) dv{y) , (3.18) 

the energy of a gas or fluid of charged particles subject to a confining potential V and 
mutual repulsion by (the two-dimensional form of) Coulomb's law. The Euler-Lagrange 
equation for the Coulomb gas energy functional reads 

V{x) -2 j \n\x - y\ du{y) + i = (a; G supp z/) , (3.19) 

where ^ is a Lagrange multiplier for the normalization constraint J dv{x) = 1. By 
differentiating once with respect to x one obtains 

V'(x) = 2P.V. f (x G supp u) . (3.20) 

J x-y 

Physically speaking, this condition means that the total force vanishes in the state of 
equilibrium inside the fluid. 

The task of determining the measure u from Eq. fl3.19p can be formulated and 
solved as a Riemann-Hilbert problem [8]. It is known that the solution is unique and 
corresponds to a minimum of the energy (hence a maximum of the integrand). 

From now on we shall simplify our work by taking the confining potential V of the 
probability measure at to be convex. This assumption ensures that the large- density 
of states is supported on a single interval: supp u = [a, b]. 

Next recall the definition ( 13. 5 p of the Cauchy transform g{z). Denoting by g± the 
two limits 

g±{x) = lim g{x ± ie) , (3.21) 

e-i-0+ 

which g{z) takes on approaching x G [a, b] from the upper or lower half of the complex 
plane, one rewrites the relation fl3.20p as 

V'{x)=g4x)+g4x) . (3.22) 

As it stands, this equation holds only for x G [a,b] C M . However, from general 
theory [9] one knows that g±{x) are the two branches of a double- valued complex- 
analytic function z {g{z), h{z)) evaluated at z = x . Thus by the principle of analytic 
continuation we have, for all z G C \ {a, b}, 

V'{z) = g{z) + h{z) . (3.23) 

Note that g{a) = h{a) = g+{a) = g-{a) and g{b) = h{h) = g+{b) = g-{b). 
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We will also make use of the integrated form of Eq. f l3.23p . For 2; 6 C \ (—00, b] let 
G{z):= I \n{z - x) dv{x) , (3.24) 

J a 

where the principal branch of the logarithm is assumed, giving G{x) G M for x G (6, 00). 
Then \et H : C \ (— 00, 6] — )■ C be defined by the equation 

V{z) =G{z) + H{z) - t . (3.25) 

We observe that G'{z) = g{z) and hence H'{z) = h{z). Moreover, by combining fl3.25p 
with the Euler-Lagrange equation fl3.19p we have 

H{b) = V{b)-G{b)+i = V{b)- f \n\b - x\ du{x) + e = G{b) . (3.26) 

J a 

By the same reasoning, H{a) = G{a). 
3.4- Asymptotics for k small 

In this subsection, we take the absolute value of the real number k to be small. In order 
to compute the large- iV asymptotics of lnQ{Nkir) for this situation, we modify the 
exact integral representation fl3.12p - fl3.16p by expressing the right-hand side of (13.160 
complex contour integral: 



where the contour Cx loops around the set of points xi, . . . , xa? . We thus obtain 

-N+l 



n{NkY]) = C2,Nk-^+^ X 
e^'^^ dz 
c. ]\^=i{.z-xi) 




\[{x,-x,f\[e-''''^-^^dxi. (3.28) 

i<j I 

Because the Coulomb gas energy is of order A^^, which is large compared to the 
perturbation due to the contour integral and here in particular the 'external electric 
field' term Nkz , we expect the density of the fluid of charges Xi, . . . , xa? to remain the 
same in the limit N ^ 00. (This will turn out to be fully correct as long as k does not 
exceed a critical value.) Thus the charges xi, . . . , xat are still expected to distribute (for 
N — )■ 00) according to the equilibrium measure v with support [a, b]. 

Taking this fact for granted, we fix a contour C encircling supp z/ = [a, b] and 
interchange the integral over {xi, . . . ,X7v} with the contour integral. Then, by taking 
the logarithm and passing to the large- limit we arrive at 

uj{k) := lim N~Hnn{NkU) 

= 7 + lim In I ^^k-N S>ik-kx)du(x) ^ ^3 29) 

with some number 7 which remains unknown for the moment, as we did not keep track 
of the overall normalization constant. Note that a;(0) = from ^2(0) = 1. 
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The integral (13.291) for — )■ cxd is computable by saddle analysis and the method 
of steepest descent. We first look for the critical points of the integrand. The condition 
for z = zq to be an extremum is 

Ja — X 

When does this equation have a solution for zq G R? To find the answer, notice that 
the function g{x) for x G M \ [a, &] is monotonically decreasing: 



X - yf 



= - I 7Zr^2 < ^ [«' ^]) • (3-31) 



We therefore have the inequality g{b) > g{oo) = > g{a) and the function 

g : R\{a,b)^[g{a),g{b)] (3.32) 

is a bijection. Thus for any k G [g{a),g{b)] there exists a unique solution 

Zo = 9-\k) (3.33) 

of the equation k = g^zo). Note that ^ 11^ \ (cl, b). 

In the following let k be fixed, with g{a) < k < g{b). To evaluate the integral (13.291) 
by steepest descent, we deform the contour C for z into the axis g~^{k) + iM parallel to 
the imaginary z-axis. Because g'{x) < for x G M \ [a, 6], the saddle at zq = g~^{k) is a 
local minimum of the integrand evaluated along the real axis, but is a local maximum 
of the integrand on the axis g~^{k) + iM . Thus the path of steepest descent leads across 
the saddle zo = g~^{k) in the direction of iM . Steepest-descent evaluation of the integral 
then yields 

uj(k) = -l + kg-\k) - f In (k g-\k) - kx)du{x) . (3.34) 

J a 

Notice that since kg~^{k) — 1 as A; — )■ , this satisfies the required normalization 
condition uj{Q) = by insertion of the additive constant 7 = — 1. For later use, we write 
our result in the equivalent form 

u{k) = -l + kg-\k) - G{g-\k)) - In . (3.35) 

Now by using -£{k z — G{z))\z=g-i(^k) = , we infer that to has derivative 

uj'{k) = g-\k) - 1/k , (3.36) 

or equivalently, 

{k-^ + u'{k)) =z. (3.37) 

k=g{z) 

Comparison with (13. 6p then shows that u'{k) = R{k) coincides with the Voiculescu 
i?-transform. Hence, by integrating, 

u{k) = / R{t)dt . (3.38) 



The reasoning above makes good sense as long as g{a) < k < g{b), so that the 
saddle zq = g~^{k) lies outside the spectrum [a, 6]. It should be mentioned that the 
same result (13.381) was established under somewhat different assumptions (see below) 
by rigorous analysis [15] using a large deviation principle. 
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3.5. Asymptotics for k large 

We turn to the complementary case of large k, meaning the ranges k < g{a) < and 
k > g{b) > 0. In this case, as we shall see, one charge dissociates from the fluid in the 
interval [a, b] in response to the strong force exerted by the external field term e^'^^. 
Consequently, the assumptions leading to the formula (13.291) are no longer met and we 
need to proceed in a different manner. 

Our modified procedure is as follows. Abandoning the contour integral formula 
(I3.27p . we insert the identities (I3.14p . (I3.16P into the eigenvalue integral representation 
(I3.12P for Q and use Sat permutation symmetry to single out one eigenvalue, say the 
first one Xi =: x . In this way we obtain the exact expression 

n{NkU) = C3,Nk-''+^ [ e^'="-^^(")7r,v-i,iv(x)rfa; , (3.39) 
7rN-i,N{x) = Z^^ I W (x, - x,f \{{x- xi) e-^^(-') dxi . (3.40) 

JRN-I 2<j<j 2<1<N 

The function 7r7v-i,Ar(x) is a polynomial of degree A^ — 1 in the variable x . By a classical 
result [H [20] it is actually the orthogonal polynomial of degree — 1 associated with the 
weight function e~^^*^^^. We find it convenient to choose the normalization constants 
Ca.AT and in such a way that HN-i.Nix) = x'^^^ + ... is monic. 

Once again, we will use saddle analysis and the method of steepest descent to 
calculate the integral (13.391) for large A^. To prepare this step, we observe that the 
orthogonal polynomial itn-i,n has a big number A^ — 1 of real zeroes concentrated in 
a small neighborhood of the interval [a,b]. The large- A^ asymptotics of the high-order 
polynomial ^n-ln divides the fc-axis into different regions. In fact, when the absolute 
value of k is large, the term e^'^^' pushes the saddle of the x-integral away from [a, b] 
and thus into the region where ttn-i,n does not oscillate but varies monotonically. The 
large- A^ asymptotic analysis of (I3.39P then is rather straightforward; see below. On the 
other hand, as k decreases below a critical value the saddle merges with the fluid [a, b], 



where h^^i^n oscillates rapidly. The integral representation (I3.39P then does not give 
a direct view of the large- A^ asymptotics. Fortunately, this case has already been dealt 
with in Section 13.41 using the alternative representation by a contour integral. 

For definiteness, from here on let k > g{b). (For k < g{a) the argument goes just 
the same.) The main contribution to the integral (13.390 then comes from large values 
of X , where the polynomial (13.401) behaves as 

71n-i,n{x) ~ e(^-l)/^°(--^)^'^(^) ~ e^^(^) . (3.41) 



By inserting this asymptotic expression into the integral (I3.39P we obtain 

uj{k) = 71 - In A; + lim A^-^ In [ ^Nkx-NV(.)+NG{x) ^3 42) 

N^oo J 

with another constant 71 . 

The condition for x = xq to be an extremum of the integrand is 

k = V'ixo)~gixo) . (3.43) 
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From fl3.23p we know that V'{z) —g{z) = h{z) is the second branch of the double- valued 
function z {g{z), h{z)). Thus the saddle-point equation is 

k = h{xo) . (3.44) 

This equation does not always have a solution for € M . However, since V is analytic 
and convex by assumption, the derivative 

h'ix) = ^ (V'(x) - gix)) = V"{x) + I' (3.45) 

for X > 6 is positive, so h{x) is monotonically increasing. Therefore ii k > h{h) = g{b), 
then k lies in the range of the function h : {b, +oo) — )■ M , and by the implicit function 
theorem the monotonicity of h guarantees the local existence of an inverse 

k ^ h-\k) = z (3.46) 

of the function z i— )■ h{z) = k . By h'{x) > for x > b the solution xq = h'^{k) of (13.441) 
is a maximum of the integrand. Hence, by applying the method of steepest descent and 
evaluating the integrand of (I3.39P at the critical point Xq = h^^{k) we obtain 

u{k)=ji + kh-\k)-Vih-\k)) + j\n du{x) 

= j^ + i + k h-\k) - H{h-\k)) -Ink , (3.47) 

where the second equality is by (13.251) . We observe that the derivative of u is still given 
by the i?-transform: 

uj'{k) = h-\k) - 1/k = R{k) . (3.48) 

By matching (I3.47P to the small-fc result (I3.34p we infer that 'ji + i = —1. (Matching 
will be justified in Section IXBl ) 

Essentially the same reasoning goes through for the opposite range k < g{a) = 
h{a) < 0, resulting in the same formula (]3.47p . Notice that the argument of the 
logarithm under the integral sign in (I3.47P remains positive in this /c-range, as there 
is a sign change in both the numerator and the denominator. 

To summarize our results, we have found that 

^ / -1 + kg-\k) - G{g-\k)) -Ink, g{a) < k < g{b) , 
^' \-l + kh-\k)-H{h-\k))-\nk, k<g{a)oig{b)<k. ^ ' 

Recalling from Section [373] the relations g{x) = h{x) and G{x) = H{x) for x G {a,b}, 
we see that the different pieces of u combine to a smooth function (Figure [I]), since each 
piece has derivative u'^k) = R{k). 



3.6. Absence of phase transition 

Related to the material of Sections 13.41 and 13.51 there exist mathematical results [15] 
by Guionnet and Maida (GM), which we now discuss briefly. Using large deviations 
techniques, GM control the large- asymptotics of the spherical integral (11. 5p under 
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Figure 1. The saddle trajectories k n- g~^{k) for g{a) < k < g{b), and k i-^ h~^{k) 
for k < g{a) or g{b) < fc, piece together at the critical points g{a), g{b) to yield a 
smooth function k i— k^-^ + R{k). 



the condition that the empirical measure (i.e., the sum of (5-functions located at the 
eigenvalues) of a sequence of diagonal N x N matrices Xn converges weakly (A^ — )■ oo) 
to a compactly supported density of states u with supp = [a,b]. In contrast to the 
present setting, GM stipulate that all eigenvalues of Xn remain strictly confined to the 
interval [a,b]. In Coulomb gas language this means that infinitely high potential walls 
are placed at the boundaries of the interval [a, b]. 

With these assumptions, GM prove that In Ixj^{k) converges to our result 
oo{k) = R{t)dt in the small-/c range g{a) < k < g{b). For k outside this range, 
however, they establish a qualitatively different behavior, separated from the small-/c 
behavior by phase transitions at the critical points k = g{a) and k = g{b). Please be 
assured that there is no contradiction with our result f l3.47p . as the underlying model 
assumptions are different. As we have seen, for large values of k the strong external 
field term e^'^^ has the effect of expelling one charge from the fluid of the other charges. 
Such a dissociation phenomenon is forbidden by the assumptions of GM. 

Still, in view of the phase transitions found by GM and the rather different integral 
representations f l3.29p and fl3.42p appearing in our treatment, the reader may wonder 
why it is that the large- asymptotics for the small-fc and large-/c regions match perfectly 
to give a smooth result for u){k). Let us now give a simple argument why this must be 
the case under our assumptions on the confining potential V. 

Consider the mapping vr : End(C^) — t- M which sends H to its projection 7r(if) = x 
to the one- dimensional subspace with projector 11 . Imagine computing the push forward 
of /xtv by this mapping. 



for large A^. In the language of quantum field theory one would call F the (large- A^ limit 
of the) generating function for the one-particle irreducible vertex functions. 



(3.50) 



Free probability meets supersymmetry 



19 



Since V is uniformly convex (i.e., its second derivative is bounded below by a 
positive constant) and analytic, so is F. Indeed, the latter results from the former 
as the fixed point of a renormalization-group type recursive process of integrating out 
variables. Now, 

uj{k) = lim N-Hn [ e^'=^"^rf/i^(i7) 

N^oo J 

= lim N'^ In [ e^'^^-^^^^^rfx = (kx - T(x)) . (3.51) 

Af-^oo J x:r'ix)=k 

Thus our function u is the Legendre transform of F. Since F(a;) is analytic and uniformly 
convex as a function of x G M, the Legendre transform uj{k) must be analytic as a 
function of A; G M, ruling out any possibility for a phase transition to occur. 

Moreover, under the stated conditions on V one can show that the /^-transform 
R{k) is an entire function of G C. Also, in the present setting there is no reason 
to think that u should have any singularities in the complex fc-plane. Hence, by the 
identity principle (i.e., uniqueness of analytic continuation) we expect that 

pk 

u{k) = / R{t)dt (3.52) 







holds for all /c G C. This is the result we propose. 



4. Fermion-fermion sector 



Our argument so far deals with the rank-one case K = (p Ip = kU . It extends to 
the higher-rank case of several replicas cpi , . . . ,ipp without difficulty. However, as we 
have seen in Section [2], the supersymmetry method requires as input the characteristic 
function Q{K) for K = Lp®Lp + %p®ip^ which includes the fermion-fermion expression 
V^®^ built from?/^ G C^®A((C^)*) and ^ G (C^)*®A(C^). The reasoning of Section 
[3] does not apply to that second summand. Therefore we now develop another scheme, 
concentrating again on the simple case of a single replica, K = ip ^ ip. 

So, in this section we make a study of Q{K) = J e^^^^^dfi^iH) for K = ip^ip. Since 
the probability measure fi^ is invariant under conjugation H i— )■ gHg~^ for g G Uat, 
the characteristic function VL{%Ij ® ip) has the property of depending only on the GLat- 
invariant {ip^ip). Although this property is not explicit from the definition of f2 , it can 
be made so by averaging ad hoc over llAr-orbits. Hence we compute 

J\Jn JVn 

as a useful preparation for the subsequent process of integrating against dfiN^H). 
We claim that the integral (14.11) has the following alternative expression: 

lN{ip®iJ]H) = N\-^ I Det{t- {ip,'^) H)e-'dt . (4.2) 

Jo 

This formula is proved in the next subsection. Having established it, we will determine 
the large- asymptotics of A^^"*^ lnQ[Nip ® ip) in Section IT2l 
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4.I. Proof of formula 

After expanding the integrand of (14.ip by the power series of the exponential function, 



Y.''^^k^9-\Hg^r , (4.3) 
^-^ ml 

m=0 

our task is to compute the integrals J^^{ipg~^, Hgip)"^dg. To do so, it is helpful to 
recall some facts from the representation theory of GL^v • 

The irreducible representations of GL^v are labeled by highest weights A which are 
sequences of integers A = (Ai, A2, . . . , Aa?) with Ai > A2 > . . . > Aat . Here we will 
need only those representations which extend holomorphically to End(C^); these are 
distinguished by Aat > 0. The character of the representation px with highest weight A 
is denoted by sx{g) = Tt pxig) and is called a Schur function. The Schur functions form 
an orthonormal system w.r.t. Haar measure dg for Uat of total mass J^^ dg = 1 : 



sx{g )sx'{g)dg = 6xx' ■ (4.4) 

To do the integrals {"ipg'^, Hgip)"^dg , our first step is to recognize the integrand 
as a special Schur function: 

( - {i^g-\ Hg^P)r = s[i'"]((^ ® i^)g-'Hg) , (4.5) 
where A = [1™] is the highest weight given by 

Ai = A2 = . . . = Am = 1 , Am+i = . . . = Aat = . (4.6) 
In the second step we use the general identity 

[ sx{Kg-^Hg)dg= [ Tr px{K)px{g'')px{H)px{g) dg 
^ Ttpx{K)Ttpx{H) ^ sxjK) sxjH) ^ 

TrpA(IdAf) SA(IdAr) 

which derives from the fact that px is a representation. 

We now establish (14. 5p . For this we start from the relation 

Det (Id;v -K) = y^ i-irs[i-]iK) , (4.8) 

»m=0 

which is a special case of what is sometimes called the dual Cauchy identity [5]. Then, 
substituting K = ip ^ ip we manipulate the left-hand side as follows: 

Det (ld;v - ^ ® V^) = e^-^ (Hiv-^®^) ^ e- E^.i — Tr m^r 

= ^+i:^^,m-^m"^ ^y"" (^,^)-. (4.9) 

' ^m=0 

By combining equations we have 

'm=0 'm=0 

Since the Schur function S[im][K) is a homogeneous polynomial of degree m in the matrix 
entries of K, the desired relation (14.51) follows on replacing ip — ipg'^Hg. 
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We are now in a position to compute our integral: 
/ i-{i^g-\Hg^)rdg 

JVn 

where we used (14. 7p . On the right-hand side we have the simphfication S[im]('?/; ^ ip) = 
{—{ip^ip))^ from (14. 5p . and the denominator is the dimension of the representation: 

SM„i(Idjv) = dim A*" (C^) = ^7^— TT . (4.12) 

' ' ^ ' m\{N-m)\ ^ ' 

Thus we obtain the following expression for 1^ : 

^ 1 /• 

1^{^®^-H) = y^ — \ ( - ^^g-\ Hgi,)rdg 
1 ^ 

■ m=0 

We finally resum this expansion. For that, we express the factorial as (A^ — m)\ = 
Jq°° t^~^e~^dt and use the dual Cauchy identity (14. 8 p in reverse to arrive at 

_ ^ (_l/t)- (H) e-'dt 

■ m=0 

POO ^ 

= N\~^ Det{t- {ip,^)) H)e'^dt . (4.14) 
Jo 

This completes the proof. 

4-2. Large-N limit 
Recall the definition 

n{Nij ®i))= I Im{N^ ®i);H) dfXNiH) . (4.15) 



As was reviewed in Section [21 in the superbosonization method one lifts Q{NiIj ® ip) 
to a holomorphic function k ^-^ Q{Nk) of a complex variable k E C . By substituting 
t — Nt in (14. 2 p and changing the order of integration, we see that 
^ j\jN+i r°° / r \ 

n{Nk) = -j^ J n Det (t - kH) dfiN{H)] e-^'dt (4.16) 

is such a function. We now tackle the task of determining its large- A^ limit. 

Stirling's formula gives limAr_5.oo A^~^ In (A^^"'""'^/A^!) = 1, and by the type of 
reasoning of Section [3] we find the large- A^ approximation 

J Det {t - kH) d^iN{H) ~ eNSp^\n{t-kx)dv{x) _ (4_^7^ 

We now take A^ — t- oo to define 

(l){k):= \im N-^\nVL{Nk) . (4.18) 

A''— >-oo 
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Note that 0(0) = 0, and by combining the above we have 

(t){k) = 1 + hm A^-Mn / q-ni+n !^\n{t~kx)du{x)^^ _ (4_^g^ 

The integral on the right-hand side is a close cousin of the integral (13.291) and, in 
fact, has essentially the same saddle point, t = k g~^{k). The overall sign in the exponent 
is reversed, but at the same time we are now integrating over t in the direction of the 
real axis, while the path of steepest descent in the case of (13.291) was along the imaginary 
direction. Thus the saddle analysis is essentially the same and need not be repeated 
here. We just quote the result: 

4>{k) = l-kg-\k)+ f In {k g~\k) - kx)du{x) . (4.20) 

J a 

This is exactly the negative of u{k) in (I3.34p . Hence we conclude (cf. (13.521) ) that 

(t){k) = - R{t)dt . (4.21) 
Jo 

Returning to the original meaning oi k = (ipjip) we state our final result as follows: 

lim N-^ In n(Ntp o ^) = - V ^ . (4.22) 

n>l 

In Section [3] we achieved a satisfactory understanding of the large- behavior of 
A^~^ lnQ{NK) for a rank-one operator K = ip ^ if. The current section does the same 
for K = ip ® ip. In both cases, our treatment carries over without difficulty to K of 
higher rank (or several replicas). 

However, the supersymmetry method calls for Q{NK) with K of the mixed type 
K = iflp + ip'^ ; cf. (12.221) . The techniques used in Sections [3] and H] are quite different, 
and at the present time we do not know how to combine them to handle the mixed 
situation. For this reason we now turn to a completely different approach. 



5. Large-A^ combinatorial theory of ^{K) 

In the present section we launch yet another attack, combining the planar graph limit 
of perturbation theory with Speicher's combinatorial approach to the free cumulants. 
An early physics paper in this direction is [28]. 



5.1. From moments to cumulants 

We begin by reviewing the usual connection between the moments and the cumulants 
of a random number, and afterwards state the analog of this connection in the setting 
of free probability theory. 
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5.1.1. Commutative case (N = 1). Although our interest is in x random matrices 
in the hmit A^ — oo , we temporarily set A^ = 1 and review the connection between 
moments and cumulants for the case of a single random variable x G M with probability 
measure /i. The moments m„ of fi, 



'rrin = / x^dfi{x) , (5.1) 
Jr 

are generated by the characteristic function, 

n{k) = / e^^'dfiix) = y"m„— , Q{0) = 1 , (5.2) 

n=o ^• 

while the cumulants c„ are generated by the logarithm of Q{k): 

OO J Yi 

lnfi(fc) = yc„— . (5.3) 

n=l 

Sometimes one wants to express the moments in terms of the cumulants, or vice versa. 
Switching between the two descriptions is an exercise in basic combinatorics and the 
Leibniz product rule of differential calculus. Let us do this exercise as a warm up for a 
more strenuous calculation to come in Section 15.51 below. 

We start by writing the n**^ moment as an ri"^ derivative at A; = : 

d 



k=o \ dk 



(5.4) 



fc=0 



By the Leibniz rule this becomes 

" \dk ^ nik) J 



(5.5) 

From this formula, by multiplying out factors using the distributive law, we generate a 
sum of terms which are in one-to-one correspondence with partitions p G n(n) of the 
set {1,2, ... ,n}. Indeed, we are to break up (or partition) n identical factors into i^iip) 
blocks of integer length / > 1 subject to Yli ^^i{p) = n . Each block of length I consists 
of one logarithmic derivative Q'{k)/Q{k) together with I — 1 derivatives d/dk acting on 
it before evaluation at A; = . Now since each block of length I contributes a cumulant 
Q = {d/dky hifl{k)\k=o , we obtain a formula expressing the moment, 

Yi (5-6) 

peU{n) 1>1 

as a sum over partitions p G n(n) in terms of the cumulants Ci, . . . ,Cn . 

5.1.2. Speicher's formula (N = oo). We return to the case oi N x N random matrices 
H with probability measure /iat and recall the definition (13.11) of the moments m„ and 
their large- A^ limits m„ = limAr_j.oo ""^n.^v • We also recall the definition (13.61) of the free 
cumulants c„ by the Voiculescu i?-transform. 
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It turns out that the large-A^ moments m„ are expressible in terms of the free 
cumulants q (with / < n) by a combinatorial formula closely analogous to (15. 6p . Due 
to Speicher [19], this formula reads 

^n= Yl n^r^'^- (5.7) 

pGNC(n) 1>1 

It differs from (15. 6 p only in that the sum over all partitions p G Il{n) has been restricted 
to the subset of non-crossing partitions p G NC(n). A partition p is called non-crossing 
if for any two pairs {^1,^2} and {ji, J2} taken from any two different blocks of p, it never 
happens that 

ii < ji <i2< 32 ■ (5.8) 

Informally speaking, this means that if we arrange the numbers 1, . . . , tt, in cyclic order 
around the boundary of a disk and connect the numbers of each block of the partition 
by lines via the interior of the disk in a 'minimal' way, then the lines associated with 
different blocks do not cross each other. See Figure [2] below for an example. 

Consider the single-block (or trivial) partition where vi{p) = . . . = Vn-i{p) = 
and z/„(p) = 1. This partition is non-crossing and contributes n«>i ^f^^^ ~ to 
Speicher's formula (15.71) for the moment m„ . Thus (15. 7p has the general structure 

mn = Cn + Pn-l(Ci, Cg, . . . , C„_i) , (5.9) 

where Pn-i is a polynomial in the free cumulants of order — 1 or less. It follows that 
the relation (15. TP can be inverted to express c„ in terms of the moments mi, . . . ,mn . 
In this sense the formula (15.70 may serve to define the free cumulants in terms of 
the large- iV moments. Speicher jl9] has proved (under suitable conditions on n^) 
that this combinatorial definition is equivalent to Voiculescu's analytic definition of 
k I— R{k) = ^Cn via inversion oi z ^ g{z) = ^ m„ z~"'~^. 



5.2. Cumulant tensor 

After these preliminaries, we will outline a combinatorial description of the large- limit 
of A^^^ \n.VL{NK) by way of expansion in powers of K (restricting the numerical part of 
K to be of finite rank). For this purpose we fix some degree n > 1 and consider the full 
tensor of cumulants 



K = 



We refer to the set of components Cj^"'jJ^ as the 'cumulant tensor' (at degree n) for 
short. Please note that our simple notation does not display the dependency which 



^ii...i„ j^g^g ^j^g random matrix dimension A^. 



For a general probability measure /ztv with no symmetries (and N > n), we expect 
all components of this tensor to be independent of each other. However, we shall 
now assume a probability distribution of the UAr-invariant form e'^^'^^^^^dH. The 
characteristic function then inherits the invariance property Q{K) = Q{g^^Kg) for 
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g eUn ■ It follows that the set of components C^^'"^" constitute an invariant tensor of 
Ujv (or by analytic extension, of GLjv = GLjv(C), the complexification of Uat): 

g G GL N ■ By a classical result of invariant theory due to Weyl [27] , one then knows that 
^ji'-'In nonzero only if the numbers {ji, . . . , jn} agree as a set with the set {ii, . . . , in}, 
i.e., if there exists an element vr of the symmetric group S„ such that ji = ij^ii) for 
/ = 1, . . . , n . Thus our cumulant tensor can be expressed as a sum over permutations: 

n 

c^:t = E^-^wn^-Mo- (5-12) 

Moreover, by its definition ( IS.lOp as an rfi^ symmetric derivative, C^^'"-^ is invariant 
w.r.t. any permutation a G S,„ of the index pairs, {iuji) ^ (iaii), ja{i)), I = 1,...,"^,. 
It follows that we may assume the coefficients 7n,Af(7i") to be conjugacy class functions, 
i.e., 7n,jv(vr) = 7n,7v(cr"^7rcr) for all a G S„ . 

Recall now from the basic theory of the symmetric group that the conjugacy class 
[vr] of an element vr G S„ is determined by its cycle structure; more precisely, by the set 
of non- negative integers z^i(vr), 1/2(7?), ... , z/„(7r) subject to Ylii ^^lij) = where I'lii^) is 
the number of cycles of vr of length /. An important role in the following will be played 
by the conjugacy class, [irr], of all irreducible cycles - by this we mean the conjugacy 
class of elements vr G S„ with z/i(7r) = . . . = z/n_i(7r) = and z/„(7r) = 1. 



5.3. Large-N hypothesis 

What has been said so far is true for all N, but now we take the large- iV limit and claim 
the following large-N hypothesis: 

hm N"-HM-) = ( ^» ''1 = ""'1 • (5.13) 

iv^oo 10 else. 

In words: A^"~"'^7n,7v(7r) goes to zero for — > cxd unless vr G S„ belongs to the conjugacy 
class of one irreducible cycle of length n . (If so, the cumulant tensor of degree n is 
determined by a single number c„ for — > 00 . We shall see that c„ is in fact the 
^th cumulant, vindicating our notation.) The hypothesis (I5.13P is present in [6], 
Theorem 4.5. To make our paper self-contained, we now offer some motivation (if not 
rigorous justification) of this hypothesis from perturbation theory. 

5.4. Perturbation theory argument 

We here assume the exponent of our probability distribution e~^^^^^^dH to be of the 
form 

Tr V(H) = -^Tr + Tr W(H) , (5.14) 
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where the parameter is small, so that the interaction (or non-quadratic part) Tr W{H) 
can be treated as a perturbation. In order to develop a perturbation expansion in the 
small parameter cr^, we single out the Gaussian (or GUE) measure 

dfiQUEiH) = e'i^^'^^ dH , (5.15) 

with Lebesgue measure dH normahzed by / dfiGVE^H) = 1. Our goal is to compute 

niNK) = J e-^^'-^(^)+^^'-^^ci/XGUE(i^) . (5.16) 

We take it for granted that n{0) = J e-^^^^^UiiQUEiH) = 1 by the choice of 
normalization constant Vr(0). 

Now, by passing to the logarithm on both sides of the equation and shifting the 
integration variable H ^ H + a'^K, we obtain 

\nn{NK) = ^Tr + In ^ e-^TrW-CH+^^K) d^i^v^iH) . (5.17) 

Next we use the trick of writing the GUE measure as the result of applying the heat 
semigroup (generated by the Laplacian A = ^ d'^/dHijdHji) to the Dirac 5-distribution 
with unit mass localized at if = : 

dfiGVEiH) = (e'^'^/2A.5)(ii) . (5.18) 

We then use partial integration to bring A^^^ \nQ{NK) into the form 

This expression serves as our starting point to develop the perturbation expansion as 
follows (we give only a brief sketch, referring to the literature [1] for greater detail). 

One expands the exponential function q-'^^^W -j^^ power series. Using 
standard graphical code, one represents each monomial Trif' in TtW{H) by an l- 
vertex. One also expands the heat operator q^^^/^^ and represents each action of the 
Laplacian A by an edge. In this way the contributions to the r.h.s. of fl5.19p are drawn 
as graphs. 

Any given graph contributes with an overall power of N^~^ where x is the Euler 
characteristic of the graph, i.e., the number of vertices minus the number of edges plus 
the number of faces. Indeed, every vertex carries a factor of A^, every edge comes with 
a factor of N~^a'^, and every face corresponds to a free summation variable, thereby 
contributing one factor of A^. (To do this counting, we first delete the vertex legs that 
are saturated by the final substitution H — > a'^K.) Thus the second term on the r.h.s. 
of fl5.19p can be organized as a sum over topological sectors: 

where u^{K) is the sum of all contributions from graphs with Euler characteristic x- 

By the linked cluster theorem, disconnected graphs cancel upon taking the 
logarithm. Now among the set of all connected graphs (with at least one substitution 
H — >■ a'^K) the Euler characteristic becomes maximal for planar graphs with the 
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topology of a disk D, where all substitutions H — )■ cr'^K are arranged around the 
boundary of D. The Euler characteristic for such a graph is x(-D) = 1; indeed, a 
triangle, say, has three vertices, three edges, and one face, so x = 3 — 3 + 1 = 1. 

The Euler characteristic of a non-planar graph is known to be smaller than the 
planar value x{D) = 1. There also exist planar graphs which have a topology different 
from that of a disk; examples are the annulus or a disk with several holes in its interior 
resulting in an 'inner' boundary. Again, the Euler characteristic of these non-disk planar 
graphs is smaller than 1. For an annulus A one has x(^) = xi^) ~ 1 = 0; since there 
is one face missing as compared with the disk D. 

The upshot of all this is that the leading contribution to f l5.20p in the large- limit 
is of order 0{N^) and is given by the sum over all planar graphs with disk topology. 
Because these graphs are constructed by inserting every substitution H — > a^K into a 
single line (circulating around the boundary of the disk-shaped graph), they all produce 
single-trace contributions TtK"^, n> 1. Perturbation theory thus leads us to expect 

oo 

lim N-Hnn(NK) =uji(K) = y"—TTK'' . (5.21) 

n=l 

The reason for writing the coefficient of Tr fT" as Cn/n will become clear shortly. 

It should be emphasized at this point that we are taking the limit N ^ oo while 
keeping the rank of (the numerical part of) K finite. Non-disk (e.g., annular) planar 
graphs make contributions of the multi-trace form 

which would not become negligible in the large- limit if we let Tr grow with A^. 

This ends our excursion into perturbation theory. Based on the result (15.211) . which 
should even be valid beyond the domain of validity of perturbation theory (if only as 
a series expansion with finite radius of convergence), it is straightforward to compute 
the large- A^ limit of the cumulant tensor (15.101) . Doing so with the help of (I5.12p . we 
immediately arrive at the large- A^ hypothesis (15.131) . Indeed, in the process of applying 
the first derivative d/ dKi^j^^^^ to Tr / n we have n factors of K to choose from and thus 
a freedom which cancels the factor 1 /n ; application of the remaining n — 1 derivatives 
d/dKi^j^^i^^ gives precisely the sum over irreducible permutations X^^gjij-r] YIi^hJ^h)- 

At the present stage, we cannot tell whether c„ indeed is the n^^ free cumulant, but 
this open question will be settled in the next subsection. Our strategy will be to show 
that the large- A^ hypothesis (I5.13P implies Speicher's formula (15. 7p . Since it is already 
known that the coefficients c„ of the latter formula are the free cumulants, the desired 
result follows. 

5.5. Retrieving Speicher's formula 

Just as in the case (15. 4p of a random number, we can generate the moments (13. ip of 
the probability measure /xat by differentiation of the characteristic function fl{K) = 
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(5.23) 



where indices are arranged in cyclic order and summed in order to manufacture the 
expectation of a single trace Trif". By the identity 



-Inn 



d 



OK, 



o e 



d 



dK„. 



+ 



dlnfl 



the expression for m^, at is cast in the form 



'v- E ( 



d 



dK 



+ 



dlnQ 



^-n^n — 1 



l-nT'Ti — l 



d 



dK 



+ 



dK 



(5.24) 



. (5.25) 



In the following we adopt the convention of assigning to d / dKi^i^_^ the number k . 

Now, by evaluating the derivatives at ii" = we get an expression which is a 
polynomial in the cumulant tensor of f l5.10p . The summands of the polynomial generated 
in this way are in one-to-one correspondence (for fixed indices and by the 

bijection d/dKi^i^^^ ^ k) with partitions p G n(n), where each block of length I of 
p corresponds to a cumulant tensor of degree /. For example, the following partition 
p e NC(8) C n(8), 

p={l,5,8}U{2,3,4}U{6,7} (5.26) 

contributes to mg, at as 



ll,...,J8 



5.5.1. Contribution from non-crossing partitions. As a first step (which will turn out 
to be the main step), we compute the contribution to f l5.25p from the non-crossing 
partitions NC(n), using the relations fl5.10p . (15.120 and the large- iV hypothesis (15.130 . 

So, let p G NC(n). To prepare for the task of counting powers of A^, we will 
associate with p a 2-complex S(p) as follows. Let D be (any) disk, and divide the 
boundary line of D into n segments numbered in counterclockwise order by 1, . . . ,n. 
These segments shall be 1-cells of the 2-complex S(p) to be constructed. Each pair 
{k,k — 1) of consecutive segments represents one partial derivative d/dKi^,i^_^, which 
we graphically depict by the boundary point (also numbered by k) between the two 
segments. These n =: points separating consecutive segments are the 0-cells of S(p). 

The partition p has not yet been used; but now, if / numbers taken from the set 
{1, . . . , ra} constitute one block of p, we draw I — 1 arcs across D to connect the members 
of that block (or rather, the numbered 0-cells assigned to them) with one another. We 
take each such arc to be another 1-cell of S(p). Note that the total number of 1-cells of 
S(p) is 

d,{p) = n + J2il-'^Mp) ■ (5.28) 

Z>1 




Figure 2. Example of a non-crossing partition p — {1,5,8} U {2,3,4} U {6,7} for 
n = 8. 0-cells belonging to the same block of the partition are connected by lines. 
These lines divide the disk into sectors (or 2-cells). 



Because the partition p is non-crossing, the hnes of the 1-cells of divide the area 
of the disk D into sectors. We take these sectors to be the 2-cells of the complex 
and denote their total number by c?2(p) . Figure [2] shows the result of this construction 
for the example given above. [Moreover, it is clear that any choice of orientation for 
the 1- and 2-cells of turns into a differential complex d) with boundary 

operator d. It follows that d) has an Euler characteristic, which can be computed 

as the alternating sum do — di{p) + d2{p)] 

We are now in a position to evaluate the contribution to fl5.25p from our fixed 
partition p G NC(n). According to the definition f lS.lOp . each block of length / of 
derivatives d/dKi^i^_^ encodes one cumulant tensor of degree / . By the identity fl5.12p 
and the hypothesis f l5.13p the large- iV leading contribution of each such block is a factor 
of N~''~^^ci times the sum over irreducible permutations J2TTe[ivv]Y[^ik,jn{k)- Note that 
the total power of from the product of these factors is 

jY-E('-i)'^i(p) — ]\fd.o-di{p) ^ (5.29) 

We still need to do the sum over indices ii,i2, . . . ,in for the product of Kronecker- 
delta symbols. Recall that by (I5.12p the set of lower indices of a nonzero component of 
the cumulant tensor must be the same as its set of upper indices. This constraint forces 
some of the indices to be equal. [For our example (I5.27P we have ii = and = ij.] 
There is a graphical meaning for these constraints: ia = ib if the two 1-cells numbered 
by a, 6 G {1, . . . , n} lie in the boundary of the same 2-cell of our complex S(p). 

Next we observe that for each block or cumulant tensor there exists just one large- 
optimal permutation vr G S; in the sum of (I5.12p - this is the shift permutation, or 
translation by one unit; it is optimal because it produces no further constraints and 
hence yields the maximal power of A^ from index sums. 

We finally calculate the index sum: each of the d2{p) 2-cells of S(p) amounts to one 
free index giving a factor of A^; there are no extra combinatorial factors, as the optimal 
permutation is unique for each block of p; hence the index sum equals N'^''^p\ It follows 




that the total contribution to ( 15.25^ from p e NC(n) is Yl cf^^'* multiphed by 

j^do-di{p)+d2{p) _ jYx(S(p)) _ (5.30) 

Now the Euler characteristic of a disk D is x{D) = 1. Therefore, since our complex 
shares with D its simplicial homology by construction, we have the relation 

= do- di{p) + d2{p) = 1 . (5.31) 

Thus the total power is A^^, which is canceled by the normalization factor in f l5.25p . 
In summary, the result of summing all contributions from non-crossing partitions 

is 

m.,^= ncr^'^ + --- (5.32) 

pGNC(n) 1=1 

where the dots signify corrections from partitions p ^ NC(n) and from the multi-trace 
terms indicated in f l5.22p . 

5.5.2. Correction terms. The corrections to (15.32^ become negligible in the limit 
— )■ oo, and we now briefly discuss why. For that, we slightly expand our graphical 
representation. Let still p G NC(n) and consider once more the associated 2-complex 
S(p). Now, however, replace the arcs for each block of p of length / by a 'blob' — 
the sum of all perturbation theory graphs for the cumulant tensor of degree / — with 
/ external lines connecting it to the I members of that block. The resulting graph is 
planar [see Figure |3] for how it looks in the case of our example f l5.27p ]. It is planar 
because the partition p was taken to be non-crossing and the blobs were chosen to be 
full disks, as opposed to disks with one or several holes in them. 

Now, inserting any correction of the multi-trace form fl5.22p amounts to changing 
the topology of a blob from disk to annulus or higher genus. Our graph then is no longer 
planar. The same goes for the replacement of p G NC(n) by p ^ NC(r2): the resulting 
graph cannot be planar, as some lines emanating from the blobs must cross. As should 
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be plausible by now, this loss of planarity results in the loss of at least one factor of A^. 
We therefore claim that all corrections vanish in the large- limit: 

n 

mn= lim m„,Ar = TTcr^^^- (5.33) 

pGNC(n) /=1 

This is Speicher's formula (15. 7p . provided that we identify the coefficients c„ - unknowns 
for us up to now - with the free cumulants of free probability theory. 

In summary, assuming that (I5.2ip holds, we have argued that the coefficients Cn 
must be the free cumulants. This concludes our perturbation theory argument in support 
of (I5.2ip . Compared to the (non-perturbative) reasoning of Sections [3] and HJ this 
argument has the advantage that it applies to the general mixed case where 

p q ^ 

K = "^^a®^a^^^h®^h- (5.34) 

a=l 6=1 

Note added. The large- hypothesis (I5.13P with Cn equal to the vl^ free cumulant, is 
essentially equivalent to Theorem 2.6 of Collins, Mingo, Sniady, and Speicher [7]. (We 
thank a referee for alerting us to that result.) 



6. Application to Disordered Scattering 

As explained in Section |2l the arrival of the superbosonization method made it possible, 
in principle, for the existing treatment of Gaussian ensembles by supersymmetry 
techniques to be extended to non-Gaussian ensembles. What was missing up to now 
was a good understanding of the large-A behavior of the characteristic function Q{NK). 
Having developed such an understanding in the present paper, we are now in a position 
to go ahead with the investigation of non-Gaussian ensembles and, in particular, of 
questions of universality. A natural candidate for a first application of the new formalism 
would be the well-known universality hypothesis for the spectral correlations of Un- 
invariant ensembles of A x A Hermitian random matrices in the large- A limit. This 
hypothesis, however, has already been discussed extensively in the literature and strong 
results have been obtained by other methods. We therefore refrain from pursuing this 
issue here and, instead, turn to a problem which so far has been inaccessible: the 
question of universality in models of stochastic scattering. 

To be concrete, we consider a standard scattering problem with M scattering 
channels coupled to an A-dimensional 'internal' space C^. We will use what is some- 
times called the Heidelberg approximation to the scattering matrix: 

S(E) = U„-2iW'^^^-j^L^W. (6.1) 

Here H is the random matrix Hamiltonian acting in the internal space C^, the scalar 
parameter E is the scattering energy, and W G Hom(C*^,C^) is a coupling operator 
with Hermitian adjoint W"^ G Hom(C^, C^^). Our interest is in the limit A — )■ cxd with 
the number M of scattering channels kept fixed. 
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To be specific, we consider the correlation function of two elements of the scattering 
matrix (correlation functions of higher order can be treated in exactly the same way): 



Cat,cd{Eu E2) = (SabiEi) - 6ab){Scd{E2) - 5,,) ) , (6.2) 



where the indices label scattering channels, and the angular brackets denote the 
expectation value w.r.t. to a probability measure /itv which we take to be UAr-invariant 
but non-Gaussian. Since we eventually want to utilize the result ( II. 4p . we require fiN 
to be of the form (11. ip with analytic and uniformly convex V. 

Let us note that the correlation function (16.21) cannot be expressed solely in terms 
of the eigenvalues of H or H^s := H — iWW'^ . One does have an alternative expression 

in terms of the Wigner-reaction matrix K{E), but there exists no transparent relation 
between the distribution of eigenvalues of K{E) and that of H. Therefore, one does not 
know how to use orthogonal polynomial techniques for the computation of (16. 2p . 

Our method of choice is the supersymmetry method as reviewed in Section |2l To 
launch it, we express the elements of the 5- matrix as derivatives of a determinant: 

Sab{E) - 5ab = -2tt|- In Det(E -H + iWXW^) . (6.3) 

aAfoa X=IAm 



For the complex conjugates Scd{E2) — 6cd in (16. 2p we do the same and thus obtain the 
following expression for the correlation function: 

Cab,cd{Ei, E2) = Z{X,Y) , 

dXbaOYdc X=Y=ldM 

f Det(Ei -H + iWXW^) Det(E2 -H- iWW^) ^ 
^) = y Det(E, -H + iWWr) Det(E, -H- iWYWf) ^^^^^^ " ^^'^^ 

While the integrand on the right-hand side differs from that of (I2.22p by the presence 
of the WW^ terms, it is easy enough to incorporate these into the formalism. We get 

Z{X, Y) = J n{K, + K2) e-T"- (^1^1+^2X2) 

^^-\{^iWW'^'fii)+\{^i,wxw'!i}i)yi{^2,WYW'^ip2)-i{i)2,ww'^i}2) 5^ 

where Kj = ipj (g) ipj + tpj ipj {j = 1,2). The integral sign stands for the Berezin 
integral over the anti-commuting variables of ijjjjipj as well as the ordinary integral 
(with Lebesgue measure) over (pj,^Pj ■ The integration domain for the latter consists of 
two copies (j = 1,2) of (C^)* x restricted to the real subspace which is given 
by (fi = —i(fi and ^2 = +iv^2 5 respectively. 

We now take the derivatives d/dXf,a and d/dYdc at X = F = Mm • We also scale 
Kj — NKj . The expression for the correlation function then becomes 

Cab,cd{E^, E2) = j n{NKl + NK2) e-^Tr(EiXi+i?,i^,) (g 

x4iVV^T^(^^-'^^)^^\V^i, W,){Wlij,){^2.W,){Wl^2). (6.7) 
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At this stage, we would like to employ the superbosonization formula f l2.27p . 
using the large- information about Q{NK) from Sections [3ll5l However, this is not 
immediately possible, as the second line of (16 .6^ is not expressed by GLAr-invariants 
owing to the presence of W and W'^ . We therefore use a trick. 



6.1. Averaging trick 

The trick is to enforce invariance by making a substitution of integration variables 
(which has unit Jacobian), 

, #i , Vj , V^i V^i ■ (6-8) 

and averaging over g with Haar measure to define the auxiliary quantity 



ab, cd ■ / 

x{i;,g-\ W,){Wlgi;,){^2g-\W,){Wlg^2)dg . (6.9) 

Our next step is to compute the unitary matrix integral D^^\^ in the large- limit. 
For this we use the following result: if A, i? G End(C^) are operators whose rank is 
kept fixed (i.e., finite) in the limit A^ — )■ oo, then 

I e-^^'^^'^'~'^dgc^J^ei-\\dN®\dN + A®B), (6.10) 

where '~' means equality in the large- A^ limit, with the right-hand side tending to the 
determinant of the Fredholm operator End(C°°) 3 X H- X + AXB . While this formula 
follows as a corollary of the relation between the integral J e"^^''*^'^^^^ ^^dg and the 
i?-transform [30], it can be obtained more directly by the observation that the matrix 
entries of g and g^^ = g"^ (under the fixed rank condition on both A and B) become 
independent complex Gaussian random variables of variance in the large-X limit. 
We now apply ( I6.10p to the present situation with A = WW'^ and B = i{Ki — K2) '■ 

^-NTr{WW''g{iKj_-iK2)g-')^g 

Uiv 

~ Det-^(Id7v ® Idjv + WW^ ® i{Ki - K2)) . (6.11) 

Then we switch from the determinant on C^(8C^ to a (super-)determinant on C^^C^'^, 
using that W'^ and W exchange with C*^, while the quadruples ipi,ip2,'ipi,'ip2 and 
(fi, (^2, ^1, ^2 exchange with the Z2-graded vector space C^'^ : 

lim f e-^T^riWWtg(iK,-iK2)g-^)^g 

= SDet-^(IdM ® Id2|2 + W^W (g) isQ) . (6.12) 
Here s = diag(l, —1, 1, —1) and Q denotes the supermatrix of GLAr-invariants 

/ i^u^i) {^1,^2) (^1,^1) (^1,^2) \ 

{^2,^1) {^2,^2) (^2,^1) {^2,1p2) 
(V^l,V22) (V^1,V'2) 
V (V^2,</'l) (V^2,¥'2) (V'2,V'l) (^2,V'2) J 



Q 



(6.13) 
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Next, to account for the post-exponential factors in the integral f l6.9p we introduce 
f{Q) = {WW ® Q) (Idiv ® Id2|2 + WW ® isQy\ (6.14) 
By a shght extension of f l6.10p to include these factors, we have 
lim d[^1^ = SDet-^(IdM ® Id2|2 + WW ® isQ) 

^4 [f{Q)ab;^i,^i f{Q)cd;^2,V2 + f iQ) cb -,^1,^2 f (Q) ad;^2 ,^1) ' (6-15) 

where the first index pair indicates the position of the matrix entry in End(C*^), and 
the second pair indicates the one in End(C^I^). In the following we write 

hm Z^ij;L=:^a.,c.(Q). (6.16) 
6.2. Superbosonization 

Since the averaging trick has converted all dependence on ip of the integrand in (16. 6 p 
into an implicit dependence through the supermatrix Q of f l6.13p . we can now apply 
the superbosonization formula (I2.27p . In the superbosonization step one forgets the 
expression (I6.13P for the supermatrix Q and treats the matrix entries of Q as the new 
integration variables; at the same time, the integrand is 'lifted' to a function of Q. 

To bring the result of superbosonization into a form suitable for the large- saddle 
analysis of Section ESI we write SDet^(g) = eNXnSHctQ ^ ^nsttIuQ^ ^-^j^ ^-^^ ^^^^^ 

motivation, we write the lift Q{NQ) of the characteristic function Q{NK) of (I5.2ip as 
fi(iVQ)=e^s^^*(«)+«(^"), $(g) = y°° (6.17) 
By sending N to infinity, we then obtain the exact result 
lim Cab,cd{z + e/N,z-e/N) 



= hm /"DQe^STr(lnQ+$(Q)-.Q)g-.STr(.Q)^^^^^^g)^ ^g^^g) 
N^oo J ' 

Here we used the fact (cf. [T7j) that the normalization constant in (I2.27P has the value 
Cp^q = 1 for the present case of p = g = 2 . The explicit expression for the invariant 
Berezin integration form DQ defined in (I2.30p becomes 

DQ = (2vr)-^ n ^^^^,^. dO~ ' ^^'^^^ 

i,j=i ^Vi,'4>j ^ipi,ipj 

We also observe that by fl2.29p the Q^^-integral ranges over (isQ)-^ = (isQ)^;^ > , 
while the Q^^-integral ranges over U2 . 

One may ask why we keep e'"^^"^^ ^^'^^ Fab, cd{Q) but neglect 0{N^) correction terms 
in Inil(NQ) which appear to be of the same order. This question will be answered at 
the end of the next subsection. 
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6.3. Saddle approximation 

In the limit N oo the dominant factor in the integral is q'^ (^'^Q+'^(Q)-^Q) _ Taking 
the first variation of the exponent gives the saddle-point equation 

oo 

Q-^ + R{Q) = z-ld2\2, R{Q) = ^CnQ''-^ . (6.20) 

n=l 

This is Voiculescu's equation ( 13. 6p . except that the role of the complex variable k is now 
taken by the supermatrix Q. 

Let 2; G M be inside the support of the large- spectral measure of the random 
matrix H . Then we know that the equation + R{q) = z for a scalar variable g G C 
has two solutions, q = g+{z) and q = g-{z), with 

^Ic g+{z) = g^{z) , 3m g+{z) = -3m g-{z) <0 . (6.21) 

The following analysis is standard [22J and we therefore give only a sketch. We first 
look for diagonal matrices Q that solve the saddle-point equation fl6.20p . The condition 
isQ^pip > selects 

Qip^ = diag {g+{z),g^{z)) . (6.22) 

(This is literally true if ^cg±{z) = 0. When the real part is nonzero, a contour 
deformation is necessary in order to reach this saddle point.) 

In the fermion-fermion sector (i.e., Q^^) there are in principle four possible choices 

of diagonal saddle for Q^^- Two of these, (++) and ( ), do not contribute for — )■ 00 

as they come with suppression factors due to fermionic Goldstone modes from 
breaking of supersymmetry. The remaining two choices, (-] — ) and ( — h), are equivalent 
in the sense that they turn out to lie on the same orbit of the symmetry group. Thus 
we select the diagonal saddle 

Qo = diag {g+{z), g^{z), g+{z), g-{z)) = ^it g+{z)ld2\2 + i3mg+{z) s . (6.23) 

The dominant part of the integrand is invariant under conjugation Q — TQT^^ by 
elements T of the Lie supergroup Ui^i|2 p]. The orbit generated by the action on Qq of 
this symmetry group is a supermanifold of saddle points 

Q = TQo{z)T-' = Dlcg+{z)ld2i2 + i3mg+{z)TsT~^ . (6.24) 

Evaluating the dominant part of the integrand along this manifold we simply get unity: 

Moreover, the integration over the Gaussian fluctuations normal to the saddle-point 
manifold also gives unity by supersymmetry. Hence we have 

lim Cab,cd{z + e/N,z-s/N) 

= j DT e-s^'-(^^««(^)^-^)F,,,,,(Tgo(^)T-') , (6.26) 

where DT denotes the Ui^i|2-invariant Berezin integration form for the saddle-point 
supermanifold. This is the very same result which is obtained for the case of the Gaussian 
Unitary Ensemble in the limit N ^ oo . In this sense the result is universal. 
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Let us summarize what are the agents of the mechanism leading to universahty. 
First of all, the dominant factor Q^^T^i^'^Q+'^iQ)-^Q) of our integral is invariant under a 
Lie supergroup Ui^i|2 • This symmetry group is determined by the type of correlation 
function under consideration and does not depend on the details of the probability 
measure fiN ■ All saddle-point supermanifolds then are orbits of Ui,i|2 . 

Second, by the principle of maximal supersymmetry the large- iV limit always selects 
the same type of orbit, Q = TQq{z)T'~^, as long as z lies in the bulk of the spectrum. [At 
the edges of the spectrum, the orbit degenerates by the vanishing of A := T'^vxg±{z).] 
The specific details of the probability measure /itv enter only via the scale factor A . 

Third, the correlation function is obtained by restricting the (non-invariant part 
of the) integrand to the orbit Q = TQq{z)T~^ and integrating along it. The scale 
factor A is needed in order to express the energy dependence of the correlation function 
in the proper units given by the mean level spacing. At the same time, the 'scattering 
observable' Fab,cd(TQQ{z)T~^) is expressed in terms of the average 5'-matrix and physical 
quantities called 'sticking probabilities' or transmission factors [22]. When this is done, 
the correlation function assumes its universal form. 

Let us finally give the reason why it was legitimate to neglect the 0{N^) corrections 
to lnQ{NQ). Any such correction is Ui^i|2-invariant and, in fact, vanishes along the 
dominant saddle-point supermanifold. Its only effect is to cause a slight perturbation 
of the scale of this supermanifold and correct the density of states by a term of order 
This effect is negligible in the large- limit. 

7. Summary and outlook 

In this paper, we employed a variety of techniques to study the characteristic function 
Q{NK) and its lift Q{NQ), which are key to a recent variant of the supersymmetry 
method, i.e., the Wegner-Efetov technique of integration over commuting and anti- 
commuting variables. What we found is that the large-A^ asymptotics of Q{NQ) for 
any unitary ensemble 'close to Gaussian' is determined by the /^-transform known from 
free probability theory. (More precisely, we made the assumption that the confining 
potential for the random matrix eigenvalues is uniformly convex and analytic, in which 
case the /^-transform is an entire function.) The task of computing correlation functions 
then reduces to a discussion centered around supermanifolds which are given as solution 
spaces of Voiculescu's equation f ll.lOp extended to the case of a supermatrix Q . 

This insight opens the door to numerous applications which in the past had been 
beyond the reach of the supersymmetry method. In the present paper, we have given 
a first application to stochastic scattering, demonstrating the universality of S'-matrix 
correlations for the case of unitary ensembles close to Gaussian. Future apphcations of 
the method will be aimed at more demanding situations with two or more cuts (i.e., 
with a density of states supported on several disjoint intervals) and double scaling limits 
at critical points. To make progress with such far-from- Gaussian problems, we first have 
to learn how to deal with singularities that develop in the /^-transform. 
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We stress that although our paper deals exclusively with unitary ensembles, the 
methods used are robust and do extend to ensembles of different symmetry type. 

Let us finish with a quick glance at a new and exciting development. In a long 
series of papers by Erdos, Ramirez, Schlein, Yau, and Tao, Vu [HI [121 [211 [131 [H] , sine- 
kernel (or GUE) universality of spectral correlations has recently been established for the 
case of Hermitian Wigner matrices, i.e., random matrices with statistically independent 
entries. Using these results as input to the present formalism, we may now address 
the wider class of random matrices given as the sum of a Wigner matrix and a unitary 
ensemble. Indeed, the characteristic function of such a random matrix is a product 
Q = ^Wigner X ^unitary In the present paper we developed a large- theory of the second 
factor, while the results of Erdos et al. give control of the first factor; more precisely, 
fiwigneriNQ) approachcs flGVE{NQ) with a rate of convergence which is sufficiently fast 
in order for sine-kernel universality to emerge. Moreover, because the sum of a Wigner 
matrix and a unitary ensemble is a sum of free random variables, the i?-transform of the 
sum is the sum of the individual i?-transforms. This offers a good prospect of obtaining 
analytical control of the more general situation. 
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